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Abstract Two-particle correlations between pions in Au+Au collisions have been 



measured at beam kinetic energies of 6, 8, and 10.8 GeV/u at the Al- 
ternating Gradient Synchrotron (AGS) over a wide range of rapidities 
using a magnetic spectrometer. The data have been analyzed in the 
Hanbury-Brown and Twiss (HBT) framework to extract source param- 
eters. The event-by-event orientation of the reaction plane has also been 
measured using a scintillator hodoscope at far forward rapidities, and 
beam vertexing detectors upstream of the target. A preliminary analy- 
sis of the dependence of the source parameters on the reaction plane is 
presented. 



Identical particle correlations have previously been used to determine 
source sizes and lifetimes of the emission regions formed in heavy-ion 
collisions [1, 2, 3, 4]. In the case of non-central collisions, such mea- 
surements usually integrate over the reaction plane of the collision and 
therefore obscure physics which depends on the relative orientation of 
the two colliding nuclei. 

In this report we discuss preliminary results from our investigations 
of the dependence of HBT source parameters on the reaction plane. 
This approach should eventually give us another perspective into the 
source dynamics. Together with other HBT dependencies of the source 
geometry and dynamics, we may thus be able to form a more complete 
picture of the nature and evolution of the emission region. 



1.1 THE HBT CORRELATION FORMALISM 



Quantum statistics require that the amplitudes for identical particles 
be added together before interpreting the square modulus as a proba- 
bility. Specifically, identical bosons, such as pions, are required to have 
symmetric two-particle amplitudes, and the HBT correlation arises from 
the cross-term in the square of the symmetrized amplitude. 

For a pair of identical bosons, emitted by an extended source /9(x) 
with four-momenta pi and p2, detected at space-time coordinates xi 
and X2, assuming that the bosonic wavefunctions can be described by 
plane waves, the symmetrized amplitude of this process is: 
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INTRODUCTION 





Which leads to the probability of detecting a pair with relative four- 
momentum q = pi — P2: 
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C2(q) ^J^J d'x,d\2 /9(xi)p(x2)|^i2p = (1-2) 

where the two-particle correlation function C2(q) is simply related to 
the Fourier transform of the source density. The normalization N(q) is 
discussed below in Section 1.4. 

The above formulation is only strictly valid for completely incoherent 
sources. To account for partial coherence effects, as well as contamina- 
tion from long-lived resonances, an empirical variable A is added to the 
definition of the correlation function as a coherence scaling parameter: 

^'^"^^ = iV(q) ^^-^^ 

1.2 HBT PARAMETERIZATIONS 

In practice, the source p{r) has usually been assumed to be Gaussian 
in configuration-space. In this case, the momentum-space distribution 
is also Gaussian: 

p(r) ~ e"^ => p(q) ~ e'^^^^ (1.4) 

For multi-dimensional HBT analyses, the relative momentum variable 
q can be expressed in terms of a variety of orthogonal components, each 
of which has model-dependent significance [5, 6]. 

In the preliminary analysis presented here, a simple 3-D Cartesian 
parameterization is chosen: {Qx, Qy, Qz), with conjugate source param- 
eters {Rx, Ry, Rz)- Rz is taken along the beam axis; Rx is lies in the 
reaction plane; and Ry is orthogonal to both. 

In this case, C2(q) has the following form: 

C2(q) = + Ae-(«-^-)'-(««-^«)'-(«^-^^)'] (1.5) 



1.3 REACTION PLANE DETERMINATION 

The reaction plane is determined using the relative orientation of two 
axes: the direction of the incoming beam particle, z, and the direction 
of the impact parameter b. In our experiment, the beam axis is defined 
by a beam vertexing detector (BVER). The EVER detector consists of 
four planes of scintillating fibers each read out by a multi-anode photo- 
multiplier tube. The fiber planes each consist of ~ 150 200 x 200 /xm^ 
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fibers, situated 5.84 m and 1.72 m upstream from the target [7]. The 
position of the projection of z onto the hodoscope can be determined 
with an accuracy of 1.5 mm at 11.4 m downstream from the target. 

Charged projectile spectator fragments are detected in a hodoscope 
(HODO), and their charge centroid calculated for each event. HODO 
consists of two orthogonal planes of 38 plastic scintillator slats with 
1 cm^ cross-sections, centered on the beam line, and situated 11.4 m 
downstream from the target. The response of individual scintillators to 
deposited charge was calibrated on a run-by-run basis and this informa- 
tion was used to find the charge-weighted centroid Q = Qxi + Qyj for 
each event, where 

The direction of the impact parameter b is then Q, defined with the ori- 
gin at the projected beam position on HODO, and the reaction plane is 
defined as the plane spanned by Q and z. x is then redefined to lie along 
Q. Implicit in this definition of the impact parameter is the assump- 
tion that the direction of proton flow - the deflection of the spectator 
fragment - is along the reaction plane. 

An estimate of the reaction plane resolution is determined by ran- 
domly dividing each event into two sub-events and looking at the {(pi — <p2) 
difference distribution for the two reaction planes calculated from each 
sub-event (see Fig. 1.1). The actual resolution for the reaction plane 
determined using the full event statistics is roughly half this value [8], 
and in this manner we obtain an estimate of 86 k, 32°. 




-150 -100 -50 50 100 150 

(f^i-tjij)'' Summed over oil Centrolities 



Figure 1,1 Distribution of the reaction plane angle difference for angles 4>i and 02 
determined by splitting each event into two sub-events. The angle resolution for the 
full event is roughly half this value {5<f) « 32°). 
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1.4 



NORMALIZATION 



The correlation function is, by definition, a normalized quantity. To 
find the proper normalization N(q), a background is generated which 
creates two-particle "events" out of single tracks from different events. 
C2(q) is then simply the ratio between real data and the event-mixed 
background. 



The shape of the measured two-particle correlation function mainly 
depends on three effects: the Coulomb repulsion between the two parti- 
cles, the two-particle resolution of the detector, and the HBT correlation 
itself. 

For the Coulomb effect, a correction fc is numerically calculated for 
a simple extended source [9, 10] of size Rq. This correction is applied 
iteratively to the background until the value of Rinv obtained from the 
one-dimensional correlation function converges to that of Rq. 

The two-particle resolution arises from the finite resolution of the 
tracking detectors in the spectrometer. Two close tracks are, at some 
point, indistinguishable from single tracks, and do not, therefore, appear 
in the measured two-particle correlation. A two-dimensional cut in rel- 
ative coordinate space, /tpr, is applied to both signal and background 
events. The separation between two tracks, projected onto the first plane 
in the spectrometer, is shown in figure 1.2, together with /tpr- 



Figure 1.2 Separation between two tracks on Tl. The /tpr cut is the sohd hne. 



1.5 



CORRECTIONS TO THE 
CORRELATION FUNCTION 




The final correlation function is fitted from a spectrum in q, which 
is generated by dividing signal events by background events. Both sig- 
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nal and background have been corrected by /tpr; the background has 
additionally been corrected by fc- 

2. RESULTS 

The reader is reminded that the following results and analysis are 
preliminary, and only contain a limited subset of the E917 data. Thus 
far, only about 10% of the data has passed through the HBT analysis. By 
December 1999, we expect at least another 50% will have been analyzed. 

The measured correlation for identical pairs of pions is shown in Fig. 
1.3, plotted as a function of Qx and Qy, where x and y are defined 
relative to the reaction plane of the event as discussed in Section 1.3. To 
improve the statistics, all of Qz has been integrated over. In addition, 
since the radii for 7r+7r+ and tt^tt" pairs are similar [10], both datasets 
were combined in the present analysis. 

Corr I 






= 2.95 +/- 


.26 




= 3.34 +/- 


.19 



Figure 1.3 \Qy\ vs. \Qx\ 
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Figure I4 \Q^\, \Qy \ < 30 MeV 



Figure 1.5 \Qy\, \Q^\ < 30 MeV 
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Central slices are shown in Figures 1.4 and 1.5 along with the pro- 
jections from the two-dimensional fit. These data have been fit witli 
the function in Eq. 1.5, and values for and Ry obtained. We find 
that = 2.95 ± 0.26 fm and Ry = 3.34 ± 0.19 fm, the latter value 
being larger than the former. The error bars are solely from the fitting 
procedure and do not include systematic errors. 

To establish a baseline for comparison, the data were also analyzed 
in the same fashion, but with x chosen along a random direction rather 
than along the reaction plane (Figs. 1.6, 1.7, 1.8). In this analysis, the 
values of R^ and Ry obtained are consistent with each other, indicating 
that the observed difference relative to the reaction plane is a real effect. 




Figure 1.6 \Qy\ vs. \Qx\ - systematic check 
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Figure 1.1 \Qx\, \Qy \ < 30 MeV 



Figure 1.8 \Qy\, |Qx| < 30 MeV 
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3. CONCLUSIONS 

In a naive geometrical model of the source, we would expect Ry to 
be larger than R^. This is consistent with the simple overlap of two 
non-centrally colliding spheres for which the overlap region resembles an 
almond in shape, with Ry > Rx- 

Other data and analyses [11, 12], however, indicate that the emission 
region is not a simple static source at these energies, with dN/dy dis- 
tribTitions and flow studies indicating both longitudinal and transverse 
expansion. The shape of the source may also be obscured by particle 
absorption by spectator matter. 

Future analyses with better statistics will allow the investigation of the 
source asymmetry with respect to the reaction plane and its dependence 
on collision centrality, pair mr, and rapidity. These, taken together with 
dN/dy distributions and flow studies, may be able to further disentangle 
expansion and absorption effects from the source shape and size, and 
thus assemble a more complete picture of the collision dynamics. This 
work is in progress. 
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